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ABSTRACT 



Context. Stellar insolation has been used as the main constraint on a planet's potential habitability. However, as more Earth-like 
planets are discovered around low-mass stars (LMSs), a re-examination of the role of tides on the habitability of exoplanets has 
begun. Those studies have yet to consider the misalignment between a planet's rotational axis and the orbital plane normal, i.e. the 
planetary obliquity. 

^ ' Aims. This paper considers the constraints on habitability arising from tidal processes due to the planet's spin orientation and rate. 

Since tidal processes are far from understood we seek to understand differences between commonly used tidal models. 
Methods. We apply two equilibrium tide theories - a constant-phase-lag model and a constant-time-lag model - to compute the 
obliquity evolution of terrestrial planets orbiting in the habitable zones around LMSs. The time for the obliquity to decrease from 
O j, an Earth-like obliquity of 23.5° to 5°, the 'tilt erosion time', is compared to the traditional insolation habitable zone (IHZ) in the 

I parameter space spanned by the semi-major axis a, the eccentricity e, and the stellar mass M s . We also compute tidal heating and 

equilibrium rotation caused by obliquity tides as further constraints on habitability. The Super-Earth G1581 d and the planet candidate 
G1581 g are studied as examples for these tidal processes. 
Results. Earth-like obliquities of terrestrial planets in the IHZ around stars with masses < 0.25 M Q are eroded in less than 0. 1 Gyr. 
Only terrestrial planets orbiting stars with masses > 0.9 M e experience tilt erosion times larger than 1 Gyr throughout the IHZ. Tilt 
erosion times for terrestrial planets in highly eccentric orbits inside the IHZ of solar-like stars can be < 10 Gyr. Terrestrial planets 
in the IHZ of stars with masses < 0.25 M Q undergo significant tidal heating due to obliquity tides, whereas in the IHZ of stars 
with masses > 0.5 M Q they require additional sources of heat to drive tectonic activity. The predictions of the two tidal models 
diverge significantly for e > 0.3. In our two-body simulations, G1581 d's obliquity is eroded to 0° and its rotation period reached its 
equilibrium state of half its orbital period in < 0.1 Gyr. Tidal surface heating on the putative G1581 g is < 150mW/m 2 as long as its 
eccentricity is smaller than 0.3. 

Conclusions. Obliquity tides modify the concept of the habitable zone. Tilt erosion of terrestrial planets orbiting LMSs should be 
included by atmospheric modelers. Tidal heating needs to be considered by geologists. 

Key words. Planets and satellites: dynamical evolution and stability - Celestial mechanics - Planetary systems - astrobiology - 
Stars: individual: G1581 - Planets and satellites: tectonics 



1 . Introduction rendering planets habitable at larger semi-major axes. Williams 

, & Kasting (1997) investigated varying obliquities for Earth and 

1.1. The role of obliquity for a planetary atmosphere concluded that a substantial part of the Earth would not be tol- 

. F* ! The obliquity of a planet, i.e. the angle ^ P between its spin axis erable for life if its obliquity were as high as 90°. Hunt (1982) 
and the orbital normal, is a crucial parameter for the possible studied the im P act of zero obliquity on the temperature distribu- 
habitability of a planet. On Earth, the orbital angular momentum tion on Earth ' s and finds a 8 lobal contraction of the inhabitable 
of the Moon stabilizes if, p at roughly 23.5° against chaotic pertur- ^ on Earth for such a case - Thus > both extremely high and 
bations from the other solar system planets (Laskar et al. 1993). ver y low obliquities pose a threat for a planet's habitability. A 
This steady tilt causes seasons and, together with the rotation d iscussion of the prospects of eukaryotic life to survive a snow- 
period of 1 d, assures a smooth temperature distribution over the bal1 Earth is S iven b Y Hoffmann & Schrag (2002). 
whole globe with maximum variations of 150 K. Measurements For iff p < 5° the habitability of a terrestrial planet might cru- 
of oxygen isotope ratios of benthicforaminifera((5 18 (9b) in deep- cially be hindered. Decreasing obliquities induce less seasonal 
sea sediments suggest that changes in the Earth's obliquity have variation of solar insolation between higher and lower latitudes, 
caused phases of glaciations on a global scale (Drysdale et al. Thus, winters get milder and summers become cooler. Given that 
2009). In general, the coupled evolution of eccentricity (e) and cool summer temperatures turn out to be more important than 
obliquity has a fundamental impact on the global climate of ter- cold winters for the emergence of continental ice sheets, smaller 
restrial planets (Dressing et al. 2010), with higher obliquities tilt angles lead to more glaciation. As a consequence, the tem- 
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perature contrast between polar and equatorial regions gets very 
strong (Spiegel et al. 2009), possibly leading to a collapse of 
the potential atmosphere (priv. comm. with Frank Selsis), which 
freezes out at the poles or evaporates at the equator. 



1.2. Variations and lockings of the obliquity 

It is generally assumed that planets form in a disk around the 
host star. It was thought that this disk is coplanar with the equa- 
torial plane of the star and the planet's spin axis is aligned with 
the orbital normal. However, observational evidence points to- 
wards a more complex formation scenario, where the stellar 
obliquity, i.e. the misalignment between the stellar rotation axis 
and the orbital plane normal, depends on stellar and planetary 
mass (Winn et al. 2010) as well as on the influence of perturbing 
bodies (Fabrycky & Tremaine 2007). Furthermore, the spin axes 
of planets can experience significant re-orientations by giant im- 
pacts during the planetary formation period. Thus, the spin axes 
of the planets themselves are not necessarily perpendicular to 
the orbital plane. The obliquity tf/ p can have any orientation with 
values < <A P < 180° (Agnor et al. 1999; Chambers 2001; 
Kokubo & Ida 2007; Miguel & Brunini 2010), where a rotation 
prograde with the orbital motion means < if/ p < 90° and 
if/p > 90° defines a retrograde planet rotation. 

In the course of the planet's lifetime, this possible spin- 
orbit misalignment can be subject to lockings or severe pertur- 
bations, e.g. by a third or more bodies inducing chaotic interac- 
tion (Laskar & Joutel 1993; Laskar et al. 1993) or Milankovic 
cycles (Milankovic 1941; Spiegel et al. 2010); the putative pres- 
ence of a moon that stabilizes the planet's obliquity as on Earth 
(Neron de Surgy & Laskar 1997); a perturber that pumps the 
planet into Cassini states (Gladman et al. 1996) or drives the 
Kozai mechanism by leverage effects (Kozai 1962; Fabrycky & 
Tremaine 2007; Migaszewski & Gozdziewski 201 1). In the case 
of G1581 d, e.g., orbital oscillations due to gravitational interac- 
tions occur with timescales of 10 3 years (Barnes et al. 2008). 

Non-spherical, inhomogeneous planets may have their spin 
state modified by gravitational interactions with other bodies in 
the planetary system, assuming the system is not coplanar. These 
mass asymmetries provide a leverage that torques the planet. The 
amplitudes and frequencies of the variations are complex func- 
tions of the orbital evolution of the system as a whole, as the 
distances between all the planets, and hence the resulting gravi- 
tational torques, are constantly changing. In many cases, such as 
on Mars, the obliquity evolution is chaotic. 

Eccentricity oscillations in the G1581 system occur on 10 3 yr 
timescales (Barnes et al. 2008). Thus, bodily obliquities will be 
significant if tidal processes between the planet and its host star 
are too weak to align the spin with the orbital plane normal on 
shorter time scales (Barnes et al. 2010b). So far, the relative in- 
clinations between most of the members of the G1581 system 
are unknown, and hence the magnitudes of the torques cannot 
be estimated at present. Only the mutual inclination of u And A c 
and d has just been measured to be nearly 30° (McArthur et al. 
2010), suggesting that large variations may be possible among 
exoplanets. Moreover, the high eccentricity of G1581 d (Mayor 
et al. 2009), likely induced by planet-planet scattering (Marzari 
& Weidenschilling 2002), suggests a significant inclination be- 
tween the orbits of planet e through f. Therefore, we should ex- 
pect that the obliquities of the G1581 planets are more likely to 
be controlled by interactions among themselves, rather than tidal 
interactions with the star. 



1.3. Tidal effects induced by planetary obliquity 

Tidal interaction between a planet and its host star may alter the 
orbital and structural characteristics of the bodies. The evolu- 
tion of the planet's semi-major axis a, eccentricity e, its rota- 
tion period P 10t ., and ifr p can put constraints on its habitability. In 
the long term, as discussed below, tidal effects reduce a planet's 
initial obliquity. We call this effect 'tilt erosion'. Given all the 
dependencies of a planet's climate on the body's obliquity, as 
described above, an investigation of tilt erosion is required. 

Tidal processes imply exchange of orbital and/or rotational 
momentum between the bodies. The total angular momentum is 
conserved but friction in the objects transforms orbital energy 
into heat, which is released in the two bodies. This effect, called 
'tidal heating', can alter the planet's geology and thus also con- 
strains habitability. On the rocky moon Io e.g., surface heating 
rates of 2 W/m 2 (Spencer et al. 2000), excited by tidal distor- 
tions from Jupiter, are related to global volcanism. It is currently 
not known if a planet's geophysical activity is controlled by the 
surface heat flux (Segatz et al. 1988; Spohn 1991; Stamenkovic 
et al. 2009; Plesa & Breuer 2009) - or vice versa. Instead, tec- 
tonic activity on a planet may be determined by the volumetric 
heating rate or the degree of lithospheric fracturing or as the case 
may be, but these properties are currently not measurable re- 
motely. While surface tidal heating rates may set an upper limit 
for a planet's habitability, they may also provide a lower limit. 
Here, we think of such planets, where tidal heating is the major 
source of inner energy that could drive geologic activity. In the 
solar system we find this constellation for various moons orbit- 
ing the giant planets, such as on Io and Europa around Jupiter, 
Enceladus orbiting Saturn, and Triton escorting Neptune. 

The role of plate tectonics for the emergence and survival 
of life has been discussed, e.g. in Kasting & Catling (2003) 
and Gaidos et al. (2005). In the terrestrial planets of the solar 
system, as well as in the Moon, the radiogenic decay of long- 
lived isotopes 40 K, 232 Th, 235 U, and 238 U provided an energy 
source that drove or drives structural convection (Spohn 1991; 
Gaidos et al. 2005). The current mean surface output of radioac- 
tive decay on Earth is « 0.04 W/m 2 while the global mean heat 
flow is « 0.086 W/m 2 (see Sect. 6 and Fig. 8 in Zahnle et al. 
2007). While an Earth-sized object obviously can maintain its 
radiogenic heat source for several Gyr, the heat flow in a Mars- 
sized planet decreases much more rapidly. Early Mars tectoni- 
cally froze when its surface heating rates, driven by radiogenic 
processes and not by tides, dropped below roughly 0.04 W/m 2 
(Francis & Wood 1982; Williams et al. 1997). Currently its sur- 
face rates due to radiogenic processes are « 0.03 W/m 2 (Spohn 
1991). Obliquity tides may deliver a source of energy sufficient 
to drive tectonic mechanisms (Williams et al. 1997) of exoplan- 
ets near or in the IHZ of LMSs. 

Besides tilt erosion, another threat for a planet's habitability 
emerges from 'tidal equilibrium rotation'. If zero eccentricity 
coincides with zero obliquity, the tidal equilibrium rotation pe- 
riod of a planet will equal its orbital period. In the reference sys- 
tem of the synchronously rotating planet, the host star is then 
static. This configuration could destabilize the planet's atmo- 
sphere, where one side is permanently heated and the other one 
dark. As we will show, synchronous rotation is only possible if 
sin(^p) = 0. 

1.4. Obliquities and the insolation habitable zone 

Previous studies (Barnes et al. 2008; Jackson et al. 2008; Barnes 
et al. 2009, 2010a; Henning et al. 2009) investigated the impact 
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of tidal heating and orbital evolution on the habitability of exo- 
planets but not the effect of tilt erosion. We investigate here the 
timescales for tilt erosion of terrestrial planets orbiting main- 
sequence stars with masses as high as 1 M Q . M dwarf stars are 
possible hosts for habitable planets (Tarter et al. 2007), they con- 
tribute by far the majority of stars in the solar neighborhood, and 
transit signatures of extrasolar planets can be distinguished most 
easily in front of small - thus low-mass - stars. We calculate tilt 
erosion times in the parameter space of e, a, and stellar mass M s . 
These regions are then compared with the traditional insolation 
habitable zone (IHZ). This zone is commonly thought of as the 
orbital distance between the planet and its host star, where stellar 
insolation is just adequate to allow the existence of liquid water 
on the planetary surface. While the first attempts of this con- 
cept (Kasting et al. 1993) were primarily based on experience 
from the solar system, recent approaches include further effects 
such as orbital eccentricity, tidal evolution, and the atmosphere 
of the planet (Selsis et al. 2007; Barnes et al. 2008; Spiegel et al. 
2009). In addition to the tidal obliquity evolution of potentially 
habitable planets, we also examine the amount of tidal heating 
induced by obliquities. 

Our models applied below describe the interaction of two 
bodies: a star and a planet. We neglect relativistic effects and 
those arising from the presence of more bodies and assume tidal 
interaction of the planet with its host star as the dominant gravi- 
tational process. Our approach does not invoke the geological or 
rheological response of the body to tidal effects. 

This article is structured as follows: After this introduction, 
we devote Sect. 2 to the two models of equilibrium tide and 
we motivate our choice of parameter values. In Sect. 3 we de- 
scribe the constraints that may arise from tidal obliquity pro- 
cesses on the habitability of terrestrial planets, i.e. tilt erosion, 
tidal heating, and equilibrium rotation. Section 4 is dedicated to 
the results, in particular to a comparison with the traditional IHZ, 
while we discuss the results in Sect. 5. In Sect. 6 we conclude. 

2. Methods 

2.1. Parametrizing tidal dissipation 
2.1.1. Equilibrium tide theories 

Most current theories for the tidal evolution of celestial bod- 
ies are derived from the assumption of equilibrium tides, orig- 
inally developed by Darwin (1879). This model assumes that 
the gravitational potential of the tide raiser can be expressed as 
the sum of Legendre polynomials P/ and it is assumed that the 
elongated equilibrium shape of the perturbed body is slightly 
misaligned with respect to the line which connects the two cen- 
ters of masses. This misalignment is due to dissipative processes 
within the deformed body and it leads to a secular evolution of 
the orbital elements of the system. The actual shape of a tidally 
deformed body is difficult to calculate directly, as it is a function 
of both the gradient of the gravitational force across the body's 
diameter, as well as the body's internal structure and composi- 
tion. Initial attempts to model tidal evolution therefore made the 
assumption that the shape of a body could be well-represented by 
a superposition of surface waves with different frequencies and 
amplitudes. This approach does a reasonable job of reproduc- 
ing the measurable tidal effects on the Earth-Moon and Jupiter- 
Io systems, but its applicability to large eccentricities (> 0.3) 
is dubious (see e.g. Ferraz-Mello et al. 2008; Greenberg 2009). 
This Tag-and-add' method consists of a Fourier decomposition 
of the tidal forcing, where small tidal phase lags siixd, or ge- 
ometrical lags 5iixi) = si(xi)/2 (MacDonald 1964), are added 



for each frequency^/. The key flaw of this theory emerges from 
the fact that the Tag-and-add' procedure of polynomials P/ is 
only reasonable for a tight range of tidal frequencies (Greenberg 
2009). 

As long as the amplitude of tidal deformation of the body 
is small, deviations from the equilibrium shape are assumed to 
be proportional to the distorting force (Love 1927). By analogy 
with the driven harmonic oscillator, the tidal dissipation function 
Q has been introduced, where 



or 1 (xi) 



tan 



{e(Xi)) 



1 ~(| -sOkT/)) ton (fifor/)) 



tan 



(e(*r/)) (i) 



(Efroimsky & Williams 2009). This parameter is a measure 
of the tidal energy dissipated in one cycle (See Appendix of 
Efroimsky & Williams (2009) for a complete discussion). 

To this point, no assumptions have been made about the fre- 
quency dependence of the lag. It is hard to constrain empiri- 
cally and depends on the structure and composition of the de- 
formed body. For the Earth, seismic data argue for Q oc with 
0.2 < a < 0.4 and yr -1 < xi ~ 10 7 Hz (Efroimsky & Lainey 
2007b, a). For gaseous planets and for stars, numerical simula- 
tions and analytical studies find a complex response spectrum, 
which is extremely sensitive to the size and depth of the con- 
vective and radiative layers (Ogilvie & Lin 2004; Goodman & 
Lackner 2009). To derive simple estimates of the effects of tidal 
dissipation on the orbital evolution, two main prescriptions for 
the lag of the equilibrium have been established over the years. 

One of these two approaches is called the 'constant phase 
lag' (CPL) model. Motivated by the slow variation of the dis- 
sipation function for the Earth, e; can be considered as con- 
stant for each tidal frequency, thus e/ oc %P (Gerstenkorn 1955; 
Goldreich & Soter 1966; Kaula 1964; Peale 1999; Ferraz-Mello 
et al. 2008). On the one hand, this procedure seems to fairly re- 
produce the rheology of the Earth. On the other hand, it results in 
discontinuities when the tidal frequency is close to zero, which 
is the case for nearly synchronized exoplanets. 

The other approach, termed 'constant time lag' (CTL) 
model, introduces the relation s/ oc ^ which is equivalent to 
a fixed time lag t between the tidal bulge and the line connect- 
ing the centers of the two bodies under consideration (Darwin 
1879; Singer 1972; Mignard 1979, 1980; Hut 1981; Eggleton 
et al. 1998). This choice coincides with the visco-elastic limit 
(Greenberg 2009). It does not produce discontinuities for van- 
ishing tidal frequencies, as they occur for the CPL model in the 
case of synchronous rotation, and it allows for a complete ana- 
lytical calculation of the tidal effects without any assumption on 
the eccentricity (Leconte et al. 2010). This is particularly rele- 
vant for the study of exoplanets in close orbits, which often turn 
out to be highly eccentric. However, the key flaw in this case is 
that coupling between the modes has been ignored. 

We use two of the most recent studies on equilibrium tides, 
one representing the family of CPL models (Ferraz-Mello et al. 
2008, FM08 in the following) and one representing the family 
of CTL models (Leconte et al. 2010, LeclO in the following), to 
simulate the tidal evolution of terrestrial planets around LMSs. 
When comparing the values of the key tidal parameters in these 
two theories, i.e. Q for the CPL model and r for the CTL model, 
we recall that both models yield similar results in the limit of 



and e — > 0. For a non-synchronously rotating object, 



the semi-diurnal tides dominate and Eq. (1) yields 
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2 \to — n\ t 



(2) 



for the leading frequency^; = 2\o) - n\. In the case of a nearly 
synchronous planet, annual tides dominate, thus xi - n an d 
Q 1 x nr (see LeclO for a detailed discussion). 

2.1.2. Response in celestial bodies 

In this paper, particular attention is drawn to the evolution of the 
spin-orbit orientation and the rotation period. As test objects we 
choose an Earth-mass planet and a terrestrial Super-Earth with 
M p = 10 MEa lt h, where M p is the mass of the planet and Meanh is 
the mass of the Earth. Terrestrial bodies of such masses may have 
a variety of compositions (Bond et al. 2010), thus their response 
to tidal processes can vary. Against this background, we choose 
a value of 0.5 for the radius of gyration r g p = I p 2 M p ^ 2 R p l of 
the terrestrial planet, where I p is the moment of inertia around 
the rotational axis. The relationship between M p and the plan- 
etary radius R p is taken as R p = (M p /ME a rth) ' 27 X ^Earth 
(Sotin et al. 2007). For both the Earth analog and the Super- 
Earth we choose a tidal dissipation value of Q = 100, owing 
to the values for Earth given by Ray et al. (2001) and Henning 
et al. (2009). Measurements of the Martian dissipation function 
QMars = 79.91 ± 0.69 (Lainey et al. 2007), as well as estimates 
for Mercury, where ^Mercury < 190, Venus, with Qvenus < 17 
(Goldreich & Soter 1966), Io, implying Q lo < 100 (Peale 
et al. 1979; Peale & Greenberg 1980; Segatz et al. 1988), and 
the Moon, with givioon = 26.5 + 1 (Dickey et al. 1994), indicate 
a similar order of magnitude for the tidal dissipation function of 
all terrestrial bodies. For Super-Earths at distances < 1 AU to 
their host star, the atmospheres are assumed to be no more mas- 
sive than 10 _1 M p (Fig. 6 in Rafikov 2006), typically much less 
massive. Although thermal tides can have an effect on a planet's 
tidal evolution (as shown for Venus by Ingersoll & Dobrovolskis 
1978), we neglect their contribution to the tidal evolution of our 
test bodies for simplicity. Over the course of the numerical in- 
tegrations we use a fixed value for Q. In real bodies, Q is a 
function of the bodies' rigidity fx, viscosity 77, and temperature 
T (Segatz et al. 1988; Fischer & Spohn 1990). A comprehensive 
tidal model would couple the orbital and structural evolution of 
the bodies since small perturbations in T can result in large vari- 
ations in Q (Segatz et al. 1988; Mardling & Lin 2002; Efroimsky 
& Lainey 2007b). To estimate the impact on our results arising 
from uncertainties in Q, we will show an example where the 
tidal dissipation function for an Earth-mass planet is varied by a 
factor of two, to values of 50 and 200, and the dissipation value 
for the Super-Earth is varied by factor of 5, to values of 20 and 
500. Henning et al. (2009) sum up estimates for Q and the Love 
number of degree 2, ki, for planets in the solar system. They find 
ki = 0.3 and Q — 50 the most reasonable choice for an Earth- 
like planet (for the relationship between Q and ko see Henning 
et al. 2009). We take Q p = lOOandfe = 0.3 (Yoder 1995) for 
the terrestrial planets, consistent with our previous studies. With 
Q' = 3/2 x Q/k 2 we thus have Q' p = 500. For the CPL model 
of FM08 we assume the dynamical Love number of the ;th body, 
£d,;, to be equal to &2,;. F° r the star we take Q s = 10 5 (Heller 
et al. 2010). 

The tidal time lag t of the Earth with respect to the Moon 
as a tide raiser has been estimated to 638 s (Neron de Surgy & 
Laskar 1997) and * 600 s (Lambeck 1977, p. 562 therein). We 
apply the more recent value of 638 s to model the tidal evolu- 
tion of the terrestrial planets with the CTL model and apply a 
fixed r s = l/(2s"ini), where «; n i is the orbital mean motion 



at the beginning of the integration. Over the course of the in- 
tegration, t p and t s are fixed. For those examples, where we 
vary Q p in the CPL model, the tidal time lag will be adapted 
tOTp = 638 s x 1 00 /Qp since roughly t p oc 1/Q p . 

Empirical values for Q p and r p are mostly found by back- 
wards integrations of the secular motions of the respective bod- 
ies, where initial conditions need to be assumed (examples for 
the derivation of Q values in the Solar System are given by 
Goldreich & Soter 1966), or by precise measurements of the 
evolution of the orbital elements. There is no general conver- 
sion between Q p and r p . Only if e = and \j/ v = 0, when 
merely a single tidal lag angle e p exists, then the approximation 
<2 P ~ l/e p ~ l/(2|n-w p |r p ) applies - as long as n-co p remains 
unchanged. Hence, a dissipation value for an Earth-like planet of 
Q p = 100 is not equivalent to a tidal time lag of 638 s, so the 
results for the tidal evolution will intrinsically differ among the 
CPL and the CTL model. However, both choices are common 
for the respective model. 



2.2. Equilibrium tide with constant phase lag (CPL) 

The reconsideration of Darwin's theory (Darwin 1879, 1880) 
by FM08 is restricted to low eccentricities and obliquities. To 
compute the orbital evolution of the star-planet system self- 
consistently, we numerically integrate a set of six coupled dif- 
ferential equations for the eccentricity e, the semi-major axis a, 
the two rotational frequencies w, (i e {s, p}, where the subscripts 
and 'p' refer to the star and the planet, respectively), and the 
two obliquities if/, given by 



de 



dt 8GM 
da a 2 



dt 4GM,M 2 



/ 49 1 \ 

— 2, A 2eo,i - yet,; + ^ + 3*5.1 O) 
- i*j > ' 

2 Z\ (4e ,i + e 2 [ - 20e ,i + -^-£ M + 



j£2,; - 3e 5 ,i] - 4 sin 2 (i/r,)[eo,/ - £8,;]j (4) 

~17 = - o,„ i n7 4g <v + A - 20e o.« + 49ei,i + s 2 ,/l+ 
dr SMirt.Rfn\ L J 

2 sin 2 (^,)[ - 2e ,i + e Sii + e 9 ,i]j (5) 

^= J^([i-6K+[i+fi]K-4- (6) 

dt 4M,r L-RfniJiX 1 7 



In these equations, Z' stands for 



* 1 



Z,' := 3G k^M JMi + Mj)-L — 
1 cP nQi 



(7) 



%i := rl ^coian/iGMj), G is Newton's gravitational constant, n 
is the orbital mean motion or orbital frequency, M, is the mass of 
the z'th body, and its mean radius. The algebraic signs of the 
tidal phase lags are given by 



Heller et al,; Tidal obliquity evolution of potentially habitable planets 



5 



where 



en,;' 










— E(2<x>; — 3n) 








= H(2a>i - n) 




'B 


£5,i 


= S(n) 




C 


eg,; 


= - 2«) 




£9,/ 


= , 


(8) 


D 



(15) 



with the sign of any physical quantity x, thus 

2(jc) = + 1 V - 1 V 0. As given by Eq. (6), the evolution of tA; 
is not an explicit function of e. However, it is implicitly coupled 
to e since dtA; = dtA;(a, <x>\), and da as well as da>\ are func- 
tions of e. Thus, for our model of coupled differential equations 
dipi = d(Ai(a(e, f), <^i(e, f))- 

Though the total angular momentum of the binary is con- 
served, tidal friction induces a conversion from kinetic and po- 
tential energy into heat, which is dissipated in the two bodies. 
The total energy that is dissipated within the perturbed body, i.e. 
the tidal heating, can be determined by summing the work done 
by tidal torques (Eqs. (48) and (49) in FM08). The change in 
orbital energy of the zth body induced by the y'th body is given 
by 



j orb.,/ 



Z' , , r 147 1 

— x [ 4e ,i + e [ - 20e ,i + — + ^e 2 ,i 



2 " 2 

- 3s 5J ] - 4 sin 2 (tA,) [s QJ - eg,,] ) , (9) 

where we define the outer braces as F 01 -b.,;. The change in rota- 
tional energy is deduced to be 

E I0t j = x ( 4e 0>i + e 2 \ - 20e ,i + 49s u + e 2 ,;l 

8 n v L J 

+ 2 sin 2 ((/r,) [ - 2s 0J + eg,; + e 9 ,i] ) , (10) 

where we define the outer braces as r rot ,,. The total energy re- 
leased inside the body then is 

EM.,i = -(E ofbJ + E roUi )>0. (11) 

If the net torque on the spin over an orbit is zero, then the body 
is said to be 'tidally locked' and the equilibrium rotation rate is 

w P qu - = «(1 +9.5e 2 ) (12) 

(Goldreich 1966; Murray & Dermott 1999; Ferraz-Mello et al. 
2008), and Eq. (11) can be written as 



pequ. 



Z' 



p -[(l + 9.5e 2 ) x r mUl r orb .,;] 



(13) 



for a planet in equilibrium rotation. Tidal heating rates in the 
planet are thus a function of e, a, and if/ p . Assuming a certain 
surface heating rate /2p qu ' = E^ p /{4jrR 2 ) is given for a planet 
in equilibrium rotation, then Eq. (13) can be solved for the cor- 
responding tilt and becomes 



|t/<p qu '| = arcsin 



-2-^ p - - (1 + 9.5e 2 ) m + C 



2® (I +9.5e 2 )+4£> 



, (14) 



For e > 0.324 Eq. (12) yields of p qa ' > 2n. Thus, the signs 
of the phase lags £n lP , si, p , s 2 , v , eg, p , and v make a discrete 
jump to + 1. Hence S — = T> and Eq. (14) has no analytic 
solution. Moreover, in this case all the terms depending on if/ p 
vanish in Eq. (11) because eo, p = eg. p and 2en,p = es iP + 
£9. p . Therefore the concept of maximum or minimum obliquity 
consistent with a given heating rate cannot apply. This situation 
is a direct result of the exclusion of higher order terms in the 
second order CPL model and is an important limitation of the 
lowe-lowi/r approximations. If more frequencies were added, this 
degeneracy would arise at higher eccentricities. 



2.3. Equilibrium tide with constant time lag (CTL) 

LeclO extended the model presented by Hut (1981) to arbitrary 
eccentricities and obliquities. For convenience and to ease com- 
parison with the CPL model, we rewrite their equations for the 
evolution of the orbital parameters as 



de 
d7 

da 
dt 

du>j 
~dt 

dt 

where 
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(20) 



£2,; is the potential Love number of degree 2 of the ;th body, and 
the extension functions in e are 
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following the nomenclature of Hut (1981). When e — all these 
extensions converge to 1. Furthermore, for r, = l/(n<2,) one 
finds Z'. = Zj. 

The tidal heating rate in the ;th body is given by 



Me) 



/3 l5 {e) /3 12 (e) 



2 — , cos(i^,) h 



1 +cos 2 (^,)-i/ 5 (e) fUi) 2 



i Me) I 



2 lf?(e) 
The equilibrium rotation rate is 

Me) 



equ. 

cl>„ = n- 



2 cos(i/r p ) 



p\e)Me) 1 + cos 2 (iA p ) 



(22) 



(23) 



(see also Levrard et al. 2007; Wisdom 2008), and therefore Eq. 
(22) can be written as 



pequ. 



P\e) 



Me)' 



fl(e) 2cos 2 ( ) A p ) 



f 5 (e) 1 + cos 2 (tA p ) 



(24) 



for a planet in equilibrium rotation. This function of ifr ? has min- 
ima at if/ p = 0° V 180° and maximum at ij/ p - 90°. The 
symmetry of tidal heating around 90° simply mirrors the fact 
that a translation between (i^ p ,Wp qu ) and {n — t/r p ,-a>p qu ') results 
in identical physical states. 

For a given surface heating rate /ip qu of a planet in equilib- 
rium rotation, Eq. (24) becomes 



3. Constraints on habitability from obliquity tides 

3.1. Tilt erosion 

We define the 'tilt erosion time' f ero , as the time required to re- 
duce an initial Earth-like obliquity of 23.5° to 5°. We numeri- 
cally integrate the two sets of coupled equations from Sects. 2.2 
and 2.3, respectively, to derive f el0 over a wide range of e, a, and 
M s . We choose two different planetary masses: an Earth-twin of 
one Earth mass and a Super-Earth of of 10 Earth masses. Planets 
are released with initial rotation periods of 1 d orbiting the IHZs 
around stars with masses up to 1 M Q . 

3.2. Tidal heating from obliquity tides 

To calculate the obliquities and \if/^ x \ of a planet we use 

the adopted minimum and maximum surface heating rates of 
0.04 W/m 2 and 2 W/m 2 , respectively. Close to the star tidal heat- 
ing will be » 2 W/m 2 even without heating from obliquity tides, 
whereas in the outer regions obliquity tides may push the rates 
above the 2 W/m 2 threshold or not - depending on the actual 
obliquity. Further outside, there will be a minimum obliquity 
necessary to yield /jp qu ' = 0.04 W/m 2 . We should bear in mind 
that these key heating rates of 0.04 W/m 2 and 2 W/m 2 are em- 
pirical examples taken from two special cases the Solar System. 
Depending on a planet's structure, composition, and age these 
thresholds may vary significantly. In general, the connection be- 
tween surface heating rates and geologic activity is still subject 
to fundamental debate. 



I* 



equ. | 



4^ 2 / lp qu y 8 15 ( e ) 



Me) 



Vl(e) 
Me) 



for4;r/? 2 / ! 77Zp > <Jy - 2jf/f 5 )//3 1 



-l/2s 



(25) 



2.4. Comparison of CPL and CTL models 

Both models converge for a limiting case, namely when both 
e — > and ij/ — > 0. If we assume that tides have circularized 
the orbit and the rotation rate of the planet has been driven close 
to its equilibrium rotation, <d p = n in the CPL model 1 , then the 
signs of the tidal phase lags become sqj = 0, = - 1, 
£2,; = 1, £5,1 = 1, £8,i = - 1, and eg; = 1. As an 
example consider the differential equation for the evolution of 
the obliquity. The other equations behave similarly. For e — > 
and iff) — > 0, Eq. (6) transforms into 



dr 2M,r 2 RhitOi 1 J ' 



(26) 



while Eq. (19) from the CTL model at first order in ip and for 
e — gives 



dr 



Z; ipi 



2M,r 2 .R 2 nu)j 



[1+6] 



(27) 



Equation (26) and Eq. (27) coincide for r, = l/(n<2,)> which is 
the value of the specific dissipation function for a quasi-circular 
pseudo-synchronous planet (see Sect. 3 in LeclO). 



1 Recall Eq. (23) and see that of^" 



also depends on if/ p in the CTL 



3.3. Tidal locking vs. equilibrium rotation 

A widely spread misapprehension is that a tidally locked body 
permanently turns one side to its host (e.g. in Neron de Surgy 
& Laskar 1997; Joshi et al. 1997; GrieGmeier et al. 2004; 
Khodachenko et al. 2007). Various other studies only include 
the impact of eccentricity on tidal locking, neglecting the contri- 
bution from obliquity (Goldreich & Soter 1966; Goldreich 1966; 
Eggleton et al. 1998; Trilling 2000; Showman 81 Guillot 2002; 
Dobbs-Dixon et al. 2004; Selsis et al. 2007; Barnes et al. 2008). 
As long as 'tidal locking' denotes only the state of dw p /df = 0, 
the actual equilibrium rotation period, as predicted by the CTL 
model of Lec 1 0, may differ from the orbital period, namely when 
e + and/or i[r p + 0. Only if 'tidal locking' depicts the re- 
cession of tidal processes in general, when e — and */^ p = 
in LeclO's model, then co p = n. As given by Eq. (23), one side 
of the planet is permanently orientated towards the star if both 
e = and i// = [ 2 ]. In this case, habitability of a planet can 
potentially be ruled out when the planet's atmosphere freezes out 
on the dark side and/or evaporates on the bright side (Joshi et al. 
1997). As long as e and if/ p are not eroded, however, the planet 
can be prevented from an o> p = n locking. In the CPL model, 
however, the equilibrium rotation state is not a function of ifr ? , 
thus 'tidal locking', denoting dw p /df = 0, indeed can occur for 
<A P * 0. 

In addition to atmospheric instabilities arising from co p = n, 
slow rotation may result in small intrinsic magnetic moments 
of the planet. This may result in little or no magnetospheric 
protection of planetary atmospheres from dense flows of coro- 
nal mass ejection plasma of the host star (Khodachenko et al. 
2007) (but see also Gomez-Perez et al. 2010; Gregory 201 1, who 



model. 



Structural inhomogeneities can also prevent rotation, as in the case 
of the Moon. 
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make counterpoints). Again, since obliquities can prevent a syn- 
chronous rotation of the planet with the orbit, its magnetic shield 
can be maintained. An example of equilibrium rotation rates is 
shown in Barnes et al. (2010a). Below, we apply the CPL and 
the CTL models to the Super-Earth G1581 d and explore its tidal 
equilibrium period as a function of e and 

3.4. Compatibility with the insolation habitable zone 

Our constraints on habitability from obliquity tides are embed- 
ded in the IHZ as presented by Barnes et al. (2008), who en- 
hanced the model of Selsis et al. (2007) to arbitrary eccentricities 
(Williams & Pollard 2002). 

The first forms of life on Earth required between 0.3 and 
1.8 Gyr to emerge (Gaidos et al. 2005). During this period, be- 
fore the so-called Great Oxidation Event * 2.5 Gyr ago (Anbar 
et al. 2007), microorganisms may have played a major role in the 
evolution of Earth's atmosphere (Kasting & Siefert 2002), prob- 
ably enriching it with CH4 (methane) and CO2 (carbon dioxide). 
On geological times it was only recently that life conquered the 
land, about 450 Myr ago in the Ordovician (Kenrick & Crane 
1997; Kenrick 2003). Therefore it is reasonable to assume that 
a planet needs to provide habitable conditions for at least 1 Gyr 
for life to imprint its spectroscopic signatures in the planet's at- 
mosphere (Seager et al. 2002; Jones & Sleep 2010) or to leave 
photometrically detectable traces on the planet's surface (Fujii 
et al. 2010). This is relevant from an observational point of view 
since the spectra of Earth-like planets will be accessible with up- 
coming space-based missions such as Darwin (Ollivier & Leger 
2006; Leger & Herbst 2007) and the Terrestrial Planet Finder 
(Kaltenegger et al. 2010a). For this reason, when our two test ob- 
jects in the IHZ show either f ero <k 1 Gyr or /ip qu » 2 W/m 2 , 
constraints on the habitability of the planet will be appreciable. 

4. Results 

4. 1 . Time scales for tilt erosion 

We plot f era for an Earth twin and a 10-Earth-mass planet orbit- 
ing a 0.1, 0.25, 0.5, and 0.75 M Q star in Fig. 1 (CPL) and Fig. 2 
(CTL). In Figs. 3 and 4 we plot f ero . for the two planets over a 
range of stellar masses and initial semi-major axes, for the two 
cases of initial eccentricities e = 0.01 and e = 0.5. We also 
include a few examples for different choices for Q p . The IHZs 
are shaded in grey. f ero , increases with increasing initial distance 
from the star as well as with decreasing initial eccentricity. In 
eccentric orbits at short orbital distances the obliquity is eroded 
most rapidly. 

For the FM08 model in Fig. 1 , we see that f ero . < 0. 1 Gyr for 
M s < 0.25 Mq for both planets in the IHZ. For M s = 0.75 M 
obliquities resist erosion for > 0.1 Gyr throughout the IHZ. 
Although we do not show the corresponding figures, we find that 
if the star has a mass > 1 M the IHZ is completely covered 
with fero. > 10 Gyr - for all the three reasonable Q values of 
both planets. 

The calculations following LeclO (Fig. 2) yield qualitatively 
similar results for small eccentricities. However, for e > 0.3 
their mismatch is of the same order of magnitude as the uncer- 
tainty in the tidal dissipation function Q (see the 0.25 M Q row). 
For these high eccentricities, the CTL model of LeclO predicts 
significantly smaller tilt erosion times, where the discrepancy is 
also stronger for lower-mass stars. 

A comparison of Figs. 1 and 2 at e = shows that, for a 
given combination of planetary and stellar mass, the contours for 



fero. = 10 7 yr meet the abscissa at almost the same semi-major 
axis. For regions closer to the star, tilt erosion times provided by 
the CTL model are shorter. Farther outside, however, the CPL 
model yields smaller values for f er o. at e — 0. On the one hand, 
these discrepancies follow from the different forms of the dif- 
ferential equations - remember that, although e — 0, the initial 
obliquity i^p = 23.5° and thus the CPL and the CTL do not con- 
verge as demonstrated in Eqs. (26) and (27) for e = A ifr ? = 0. 
On the other hand, the incompatible values for Q p and r p (see 
Sect. 2.1.2) lead to different evolutionary speeds in the two mod- 
els. 

The projection of f el0 onto the M s -a plane based on the 
FM08 model is shown in Fig. 3. For the Earth twin in the 
e = 0.01 case (upper left panel) we find that the complete IHZ 
is covered by t ao . > 1 Gyr as long as M s > 0.8 M Q , where 
the inner border of the IHZ is at roughly 0.4 AU. For an eccen- 
tricity of 0.5 (lower left panel) the Earth-sized planet requires 
an M s > 0.85 M Q host star, thus a minimum semi-major axis 
of 0.5 AU, for tilt erosion times larger than 1 Gyr throughout the 
IHZ. The larger moment of inertia of the 10ME art h Super-Earth 
(panels in the right column) assures that the initial obliquity of 
23.5° is washed out on timescales > 1 Gyr even for slightly 
tighter orbits in the IHZ and somewhat lower-mass stars at both 
eccentricities. Comparing the upper and lower panel rows, cor- 
responding to e — 0.01 at the top and e = 0.5 at the bottom, 
one sees that higher eccentricities tend to erode the spin faster 
by a factor of 3. 

The higher order terms of e in the equations of LeclO pro- 
duce a large discrepancy in f ero . between e = 0.01 (upper 
panel) and e — 0.5 (lower panel) in Fig. 4. While the re- 
sults for e = 0.01 almost coincide with those of the FM08 
model, LeclO predicts much more rapid tilt erosion for e = 0.5. 
The magnification of tidal effects at high eccentricities is due 
to the plummet of the mean orbital distance in eccentric orbits 
and is enhanced by the steep dependency of tidal effects on a 
(Wisdom 2008; Leconte et al. 2010). This behavior at high ec- 
centricities is due to the slow convergence of elliptical expansion 
series and has been studied in detail by Cottereau et al. (2010). 
For the lower left panel, where e = 0.5 for the Earth-sized 
planet, an M s > 0.9 M host is necessary for f ero . > 1 Gyr over 
the whole IHZ, with the inner border at about 0.55 AU. Again, 
the Super-Earth is slightly more resistive to tilt erosion. A com- 
pelling outcome of this visualization is that the LeclO model 
predicts f ero < 10 Gyr for an Earth-mass planet at the inner 
edge of the IHZ of a 1 M star if the initial eccentricity is larger 
than 0.5. 

4.2. Tidal heating from obliquity tides 

In Fig. 5 we plot ifs^^ and as a function of a and e for 
an Earth-mass planet orbiting stars with four different masses, 
as above. Results for the 1 OA/earth terrestrial planet are not 
shown because they are very similar. The left column shows 
the results for the CPL model and the right column the CTL 
model. The white space in the left panels is due to the am- 
biguity of Eq. (14) as explained at the end of Sect. 2.2. Red 
zones indicate that even for the case of if/^ < 1 °, tidal sur- 
face heating rates are /ip qu ' > 2 W/m 2 . Starting from the left 
in a certain panel, the contours illustrate a maximum obliq- 
uity of 1°, 20°, 30°, 40°, 50°, 60°, 70°, 80°, and 89°, all of which 
produce /jp qu ' = 2 W/m 2 at their respective localization in 
the a-e plane. Green in Fig. 5 depicts a region where neither 
any obliquity produces tidal surface heating rates > 2 W/m 2 
nor any minimum obliquity is required to yield tidal surface 
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heating rates > 0.04 W/m 2 . The 9-tuple of contour lines for 
ifr^ refers to obliquities of 1°, 20°, 80°, and 89° providing 
/jp qu ' = 0.04 W/m 2 at their respective localization. Smaller 
obliquities lead to less tidal heating. In the blue zone finally, not 
even a planet with a rotational axis perpendicular to the orbital 
plane normal, e.g. ip p = 90°, yields /ip qu ' ss 0.04 W/m 2 . 

Both models predict that Earth-like planets in the inner IHZ 
of stars with masses as 0.1 M Q are subject to intense tidal heat- 
ing, while tidal heating in the outer zone of the IHZ is more mod- 
erate for the CPL model. With increasing stellar mass the IHZ is 
shifted to wider orbits for both models, whereas the green stripe 
of moderate tidal heating is located at roughly the same position 
in all the 8 panels of Fig. 5. This is due to the weak sensitivity 
of i/<p qu on M s compared to the sensitivity on a. The IHZ of a 
0.25 M star nicely covers the zone with moderate tidal heating. 
In this overlap, adequate stellar insolation meets tolerable tidal 
heating rates. Only highly eccentric orbits at the inner IHZ are 
subject to extreme tidal heating, as shown in the right column 
for the CTL model. Terrestrial planets in the IHZ of 0.5 M stars 
do not undergo intense tidal heating. Only at the inner border of 
the IHZ /Zp qu ' ~ 1 W/m 2 for high obliquities, while at the outer 
regions tidal heating rates are of order 10mW/m 2 and smaller 
even for tf/ p = 90°. For terrestrial planets in the IHZ of stars 
with masses > 0.75 M , both models indicate that tidal heating 
has negligible impact on the planet's energy budget and that it 
does not induce constraints on the planet's habitability. 

4.3. GI581 of and g as examples 

G1581 d is an M p * 5.6 AfEanh Super-Earth, grazing the outer 
rim of the IHZ of its M s =a 0.29 M Q host star (Bonfils et al. 
2005b; Selsis et al. 2007; von Bloh et al. 2007; Beust et al. 
2008; Barnes et al. 2009; Mayor et al. 2009; Wordsworth et al. 
2010; Vogt et al. 2010). The apoastron is situated outside the IHZ 
whereas the periastron is located inside. While the solar con- 
stant, i.e. the solar energy flux per time and area on the Earth's 
surface, is about 1000 W/m 2 when the Sun is in the zenith 3 , the 
stellar incident flux on G158 1 d averaged over an orbit (Williams 
& Pollard 2002) is 

/ = Lg ' 581 * 432 W/m 2 , (28) 

4tto 2 Vl - e 2 

with the luminosity of the host star Lgi58i = 0.013 L (Bonfils 
et al. 2005a; Mayor et al. 2009). If the absorbed stellar flux per 
unit area is larger than ss 300 W/m 2 , runaway greenhouse ef- 
fects may turn a terrestrial planet uninhabitable (Zahnle et al. 
2007; Selsis et al. 2007, and references therein). Thus, the max- 
imum bond albedo for G1581 d compatible with this habitability 
criterion is 300/432 as 0.7. A recent study on atmospheric con- 
straints on the habitability of G1581 d is given by von Paris et al. 
(2010). 

If confirmed, the recently detected planet candidate G1581 g 
(Vogt et al. 2010; Anglada-Escude 2010) is located more suit- 
ably in the IHZ of G1581, but it does not show a significant ec- 
centricity. Since eccentricity has a crucial impact on a planet's 
orbital evolution, we choose G1581 d as an example to study the 
change of its spin orientation and period as well as its current 
tidal heating. An analysis of tidal heating in G1581 g is given at 
the end if this section. Since we take the parameters for G1581 d 
from Mayor et al. (2009) and the parameters for planet g from 

3 Outside of the Earth's atmosphere, i.e. at the distance of 1 AU from 
the Sun, the solar constant is about 1400 W/m 2 . 



Vogt et al. (2010), our model for the two tidal evolutions need 
not necessarily be consistent. 

For the tidal evolution of G158 1 d, a comparison between the 
FM08 and the LeclO models is applied where it is reasonable. 
As initial obliquities we take 23.5° for the planet and 0° for the 
star, the initial rotation periods are 1 d for the planet and 94.2 d 
for the star (Vogt et al. 2010). The remaining physical quantities 
of the host star and the two planets are assumed as explained in 
Sect. 2.1.2. For both models we find insignificant decreases in 
e during the lifetime of the system, which has been estimated 
to > 2 Gyr (Bonfils et al. 2005b) and « 4.3 Gyr (Vogt et al. 
2010). In the long term, FM08's model predicts a negligible 
increase of a, while LeclO's model yields an insignificant de- 
crease. As shown in the left panel of Fig. 6, both models predict 
a steady increase of the rotation period of G158 1 d. For the model 
of LeclO, the numeric solution for the equilibrium rotation pe- 
riod of P e qu. ~ 35 d can also be calculated analytically with 
Eq. (23), using the observed eccentricity e = 0.38 (± 0.09) and 
assuming that the spin axis is co-aligned with the orbit normal, 
kjj p = 0°. In the right panel of Fig. 6, both models predict d's 
obliquity has been eroded over the course of the system's evolu- 
tion, assuming no perturbations from other bodies. At an age of 
< 0.1 Gyr, tidal processes would have pushed an Earth-like ini- 
tial obliquity to zero and the rotation is most likely caught in its 
equilibrium state. In both models the obliquity of the planet first 
increases to « 51° before it is finally eroded. The reason for the 
change in the algebraic sign of di/r p /d? is found, when one sets 
d^p/df > 0. Considering that £ p «c 1 for G1581 d, dtA p /df > 
in Eq. (19) is equivalent to 

w p 2f 2 (e) 

— > -— (29) 

n ~ f 5 (e)/3 3 (e) cos^p) ' 

For iff p = 51° we find oj p /n > 6.13. After about 25.3 Myr 
in our simulations with the model of LeclO, the rotation has 
slowed down to u> p /n < 6.13 (see left panel in Fig. 6), thus 
di/ip/df < 0. For the rest of the evolution dt/^p/df < since the 
orbital period keeps increasing until it reaches the equilibrium 
state. Comparing both panels of Fig. 6, we find that the time 
scale for the relaxation of the obliquity is longer than the time 
scale for the restoration of the equilibrium rotation period. This 
behavior has recently been explained by Correia (2009). 

As given by Eq. (23) from the CTL theory and as explained 
in Sect. 3.3, the equilibrium rotation period is also a function of 
e and kjj p . We calculate the equilibrium rotation period P equ as a 
function of the obliquity ip p for a set of eccentricities (Fig. 7). 
The most likely value of P equ ss 35 d, corresponding to 
e = 0.38 and i/r p < 40°, can be inferred from the gray line in 
Fig. 7 for e = 0.4. Since the observed orbital period of G1581 d 
is about 68 d, Fig. 7 shows that G1581 d does not permanently 
turn one hemisphere towards its host star, except for the cases of 
ijj p ss 74° V 106°. And even then the obliquity is so large that 
most of the surface cycles in and out of starlight over an orbit. 

We have done the same calculations of the tidal evolution 
for G1581 g, assuming a slight orbital eccentricity of 0.01 and a 
mass of 3 . 1 ME- d nh, corresponding to its projected minimum mass 
(Vogt et al. 2010). In brief, we find an evolutionary scenario sim- 
ilar to that of planet d but on a considerably shorter timescale. 
After as 20 Myr both tidal models predict a locking in equilib- 
rium rotation at almost exactly the orbital period of ~ 36.56 d. 
The evolution of ijj p for planet g looks almost the same as for 
planet d, except that the obliquity is eroded about a factor of 5 
faster in FM08's model and about a factor of 3 faster in the model 
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of Lec 10, due to the smaller planetary mass and semi-major axis 
of planet g. 

In Fig. 8 we present the tidal surface heating rates, includ- 
ing obliquity heating, of G1581d and g as a function of their 
putative obliquity. Since these planets encounter considerable 
interactions with the remaining four planets of the system, all 
orbits might be significantly tilted against each other (see Sect. 
1.2). Thus, the obliquities \jj p of the planets might be significant. 
Between spin-orbit alignment (ifr p = 0°) and a planetary spin 
axis perpendicular to the orbital plane (i// p = 90°) tidal sur- 
face heating rates on G1581 d (left panel in Fig. 8) range from 
roughly 2 to about 15mW/m 2 , given that e = 0.38 (+ 0.09). 
As a comparison, we show the tidal surface heating of the po- 
tentially habitable planet G1581 g in the right panel of Fig. 8. 
Measurements by Vogt et al. (2010) indicate that the eccentric- 
ity of G1581 g is small if not zero. As long as its eccentricity is 
< 0.3 we find tidal heating on G1581 g to be < 150mW/m 2 . 
Only if the planet's eccentricity has been > 0.5 in the history 
of the system, high obliquities could have produced tidal sur- 
face heating rates of « 1 W/m 2 (upper half in the right panel of 
Fig. 8). 

The eccentricities of the planets around G1581 could be sub- 
ject to fluctuations. In their history, their values might have been 
significantly larger than their current values. For reasonable ec- 
centricities < 0.7 of G1581d, we find that tidal surface heat- 
ing rates on this planet have always been smaller than 1 W/m 2 if 
the planet's semi-major axes has never been significantly smaller 
than its current value. On planet g tidal heating will be negligible 
if its eccentricity turns out to be < 0.3. 

5. Discussion 

Our results for tilt erosion, tidal heating, and equilibrium rota- 
tion picture scenarios for atmospheric scientists and geophysi- 
cists. For low planetary obliquities, atmospheres of planets at the 
outer border of the IHZ can be subject to freeze out at the planet's 
poles (priv. comm. with Franck Selsis). Moreover, tfp — > 
weakens the seasonal temperature contrast at a given location 
on the planetary surface but strengthens the temperature contrast 
as a function of latitude (see Sect. 1.1). Given the mass and the 
age of a star, as well as the orbital parameters of the planet, our 
estimates of the tilt erosion time scales indicate whether a sig- 
nificant obliquity of the planet is likely. Our choice of 5° as the 
final obliquity for the computation of f ero . is somewhat arbitrary 
but the time for i[r p to decline from 5° to 0.5° and even to 0.05° 
is usually much smaller than f e ro. as we use it (see right panel in 
Fig. 6). In any case, if low obliquities constrain the habitability 
of exoplanets, then the respective effects are likely to occur for 
tffp of the order of a few degrees and far before it has reached 0. 

The spin angular momentum history of a planet, i.e. the evo- 
lution of obliquity and rotation rate, is a crucial feature of plan- 
etary climate. We encourage climate modelers to consider the 
diversity of atmospheres that are possible within the constraints 
outlined above. 

Our exploration of tidal heating will be relevant for struc- 
tural evolution models of terrestrial planets. Moreover, the tidal 
response is a function of the internal structure, hence a natu- 
ral feedback mechanism exists between the two. The structure 
and cooling of Super-Earths have recently been explored in the 
literature (e.g. by Sleep 2000; Valencia et al. 2007; O'Neill & 
Lenardic 2007; Kite et al. 2009; Stamenkovic et al. 2009; Plesa 
& Breuer 2009), but a consensus has yet to emerge. On the one 
hand, the heat flow out of the interior should scale as the ra- 
tio of volume to the surface area, i.e. oc R h This relationship 



just assumes that larger planets have more radiogenic isotopes. 
Therefore large terrestrial planets seem more likely to have mo- 
bile lids than the Earth. On the other hand, larger planets have 
stronger gravitational fields and hence planets may fuse together 
due to the higher overburden pressure. Therefore large terrestrial 
planets seem more likely to have stagnant lids than the Earth. As 
if the story was not complicated enough already, none of those 
studies considered any potential tidal heating. Hence the issue 
for planets in the IHZ of LMSs is further complicated, but at least 
appears ripe for observational constraints: As recently shown by 
Kaltenegger et al. (2010b), acute volcanism on Earth-like plan- 
ets will be detectable with the James Webb Space Telescope. Our 
calculations suggest that such worlds orbiting < 0.25 M Q stars 
orbiting closer than 0.1 AU will present promising targets. 

Tidal heating rates, as given by the models of Lec 10 and 
FM08, are a function of e and if/ p , amongst others. For G1581 d 
wefind2mW/m 2 < h^^Jjf/p) < 15 mW/m 2 for its current 
eccentricity. This is in line with the value of 0.01 W/m 2 esti- 
mated by Barnes et al. (2009), who neglected the planet's obliq- 
uity. We also explore /j^g^ as a function of the planet's eccen- 
tricity and find that tidal heating is unlikely to have constrained 
the habitability of the planet since h^ si d < 0.1 W/m 2 as long 
as e < 0.7. Tidal heating on the exoplanet candidate G1581 g, 
which would orbit its host star in the IHZ, is < 150mW/nr as 
long as its eccentricity is smaller than 0.3. 

Selsis et al. (2007) approximated the rotational relaxation 
time of the planet to be lOMyr, concluding that the planet has 
one permanent day side and night side. Indeed, with the revised 
value of e = 0.38 (+ 0.09) (Mayor et al. 2009) we find that 
the tidal equilibrium rotation occurs after approximately 80Myr 
(Lec 10) - 90Myr (FM). However, the equilibrium rotation pe- 
riod of the planet, as calculated with LeclO's ansatz, turns out 
tobeP^5 gld * 35 d * P orb ./2. Correia et al. (2008) have 
approached the issue of equilibrium rotation of terrestrial plan- 
ets for low eccentricities and obliquities. For G1581 d they found 
cf qu - In * 125oP equ - * 54 4 d 

If the eccentricity of G1581 g turns out to be and if this 
state was stable at least over the recent 1 Myr, then tides have 
driven this planet into an equilibrium rotation state, where both 
\j/ p = and e — 0. Due to the smaller planetary mass and 
semi-major axis, as compared to planet d, G1581 g is driven into 
its equilibrium state faster by a factor between 3 (FM08) and 5 
Lec 10 in t ero .. 

The semantics of the terms 'tidal locking' and 'synchronous 
rotation' carry the risk to cause confusion. Wittgenstein (1953) 
has pointed out that the meaning of a word is given by its us- 
age in the language (see §30, 43, and 432 therein). However, 
'tidal locking' and 'synchronous rotation' are used ambiguously 
in the literature (see Sect. 3.3), thus their meanings remain dif- 
fuse. While some authors simply refer to 'tidal locking' as the 
state where a planet permanently turns one hemisphere towards 
its host star, others use it in a more general sense of tidal equilib- 
rium rotation. Furthermore, the term 'pseudo-synchronous ro- 
tation' appears, depicting dw p /d? = while u> p + n. We 
caution that in the CTL model of Lec 10 tidal processes drive 
the planetary rotation period towards the orbital period only if 
e — and ifr p = 0. If a planet with e + or i/f p + is 
said to be tidally locked, then it is clear that tidal locking does 
not depict the state of the planet turning one hemisphere perma- 
nently to its star. However, for certain combinations of e + 
and tf/ p + the equilibrium rotation period can yet be equal 
to the orbital period (see Fig. 7 for G1581 d, where <d p = n if 
\jj p a; 74° V 106°). And in general, the rotation period of the 
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planet will not be synchronized with respect to its orbit, if 'syn- 
chronous' here means 'equal' or 'identical'. In general, we call 
dwp/df = simply the state of 'rotational equilibrium' . When 
rotational equilibrium meets ai p + «, we call this constellation 
'pseudo-synchronous'. 'Tidal locking' is used when there is no 
more transfer of angular momentum over the course of one orbit, 
i.e. e = 0, iff p = 0, and a> p = n. Strictly speaking, a body with 
e + or^p + cannot be tidally locked because e and if/ p 
evolve with time. 

Uncertainties in our calculations arise from the assumption 
of a constant tidal dissipation in the planets, parametrized by the 
tidal dissipation factor Q or t. For one thing, merging all the 
geophysical effects such as composition, viscosity, temperature 
distribution, and pressure in one parameter unquestionably is an 
oversimplification. In addition, it is unknown how Q depends on 
the respective tidal frequency. Although on Earth Q is constant 
over a a wide range of frequencies, the behavior will be differ- 
ent for different objects. The Q value of terrestrial exoplanets 
remains to be estimated, either by observational constraints or 
theory. For another thing, whatever Q value will turn out to de- 
scribe terrestrial planets most reasonably, it will be a function of 
time due to the structural and rheological evolution of the planet. 

Our results need to be considered against the background of 
further gravitational processes, such as the entourage of mas- 
sive moons, the Kozai mechanism, Cassini states, planet-planet 
scattering, and inhomogeneities in the structure of the invoked 
bodies. These states and processes will add to - and eventually 
counteract - the tidal evolution. In those cases, significant obliq- 
uities could be maintained on time scales much larger or shorter 
than we present here. 

6. Conclusions 

Tidal processes raise severe constraints on the habitability of ter- 
restrial planets in the IHZ of LMSs. First, tidal erosion of the 
obliquities, i.e. 'tilt erosion', of such planets occurs much faster 
than life emerged on Earth. Planets with masses < 1 OA/garth or- 
biting < 0.25 M Q stars in the IHZ lose initial Earth-like spins of 
23.5° within less than 0.1 Gyr, typically much faster. In order to 
cover the complete IHZ with t ao , > 1 Gyr at high eccentrici- 
ties, stellar masses > 0.9 M are required. For terrestrial planets 
formed in highly eccentric orbits inside the IHZ, tilt erosion can 
occur within 10 Gyr up to stellar masses of 1 M . 

Second, tidal heating of Earth-like planets in the IHZ of stars 
with < 0.25 M Q is significant. Depending on the body's obliq- 
uity, orbital semi-major axis, eccentricity, and its structural com- 
position, a terrestrial planet in the IHZ of a < 0.25 M Q host star 
can undergo intense tidal heating with tidal surface heating rates 
» 1 W/m 2 . We find tidal surface heating rates of 2mW/m 2 
to 15mW/m 2 , depending on if/ p , for the Super-Earth candidate 
G1581 d and < 150mW/m 2 forG1581 g. Simulations of coupled 
orbital-structural evolution of Earth-like planets are necessary to 
further explore the effect of tidal heating on habitability. 

Third, the tidal equilibrium rotation period provided by 
LeclO is a function of the planet's eccentricity and its obliquity. 
In general, as long as the planet maintains a significant obliq- 
uity or eccentricity, its rotational equilibrium period will not 
match the orbital period. As an example, we deduce the equilib- 
rium rotation period of the extrasolar Super-Earth G1581 d to be 
around 35 d. With an orbital period of roughly 68 d, this means 
^Gis'sid ~ A>rt>./2. The model of FM08 neglects the impact of 
obliquity on the equilibrium rotation state. 

For low eccentricities (e < 0.3) and low obliquities (i/r p < 
a few degrees), the CPL model of FM08 and the CTL model of 



LeclO mathematically converge. In our calculations of f el0 (Sect. 
4) they provide similar results in the low-eccentricity regime. 
Since the model of LeclO considers higher terms in e, tidal pro- 
cesses in the moderate and high eccentricity regime occur on 
significantly shorter time scales than predicted by FM08. Hence, 
tidal heating computed with LeclO is also more intense. 

The possible detection of exomoons via transit photometry 
(Sartoretti & Schneider 1999), the Rossiter-McLaughlin effect 
(Simon et al. 2010), planet-moon eclipses (Cabrera & Schneider 
2007), and transit timing and duration variations (Szabo et al. 
2006; Kipping 2009) as well as the empirical constraint of the 
oblateness of a transiting planet (Carter & Winn 2010a,b), are 
first steps towards the measurement of the obliquity of exoplan- 
ets. All these observational effects are only accessible for transit- 
ing planets. Since obliquity determines atmospheric conditions 
as well as the amount of tidal heating and the equilibrium rota- 
tion period, it will be indispensable to verify tf/ p for planets in 
the IHZ around their host stars to assess their habitability. Thus, 
transiting planets are the most promising targets for a compre- 
hensive appraisal of an extrasolar habitat. 
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Fig. 1. Tilt erosion times for FM08. The IHZ is shaded in gray, contours of constant f ero are labeled in units of log(f ero ./yr)- Error 
estimates for Q p are shown for the 0.25 M Q star. 
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Fig. 2. Tilt erosion times for LeclO. The IHZ is shaded in gray, contours of constant f ero are labeled in units of log(f elo ./yr)- Error 
estimates for Q p are shown for the 0.25 M Q star. 
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Fig. 3. Tilt erosion times for FM08. The IHZ is shaded in gray and the contours of constant tilt erosion times are labeled in units of 
log(f ero ./yr). 
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Fig. 4. Tilt erosion times for LeclO. The IHZ is shaded in gray and the contours of constant tilt erosion times are labeled in units of 

log(fero./yr). 
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Fig. 5. Obliquity thresholds i// mdx and if/ m m as explained in the text. For each 9-tuple of contours the lines indicate 
1°, 20°, 30°, 70°, 80°, and 89°. The IHZ is indicated with dashed lines. In the red zone, h e p qu - > 2 W/m 2 for any obliquity. 
For e > 0.327 the CPL model breaks down, so this region is left blank. 
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Fig. 6. Left: Evolution of the putative rotation period for G158 1 d for both models, FM08 and LeclO. The initial period was taken to 
be analog to the Earth's current day. The filled circle at » 25.3 Myr on the LeclO track indicates w p /« < 6.13, thus from then on 
di^p/df < 0, as given by Eq. (29). Right: Evolution of a putative initial Earth-analog obliquity for G1581 d for both models, FM08 
and LeclO. 



100 




Fig. 7. Equilibrium rotation period of G1581 d as a function of obliquity for various eccentricities, following the model of LeclO. 
The gray line corresponds to e — 0.4, close to the observed eccentricity of e — 0.38 (+ 0.09) (Mayor et al. 2009). The observed 
orbital period is indicated with a dotted line. 
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Fig. 8. Tidal surface heating rates, including eccentricity and obliquity heating, on G1581 d (left panel) and G1581 g (right panel) 
following the model of LeclO as a projection on the parameter plane spanned by obliquity tfr p and eccentricity e. Equilibrium rotation 
of the planets is assumed. The observed eccentricity of 0.38 (± 0.09) for G1581 d (Mayor et al. 2009) is indicated with a dashed line, 
while the eccentricity of G1581 g is compatible with 0. Depending on obliquity, tidal surface heating rates on G1581 d vary between 
roughly 2 and 15 mW/m 2 , while on planet candidate g they are < 150 mW/m 2 as long as e < 0.3. 



